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Computing  the  Reliability  of  k  out  of  n  systems 

by 

Philip  J.  Boland  and  Frank  Proschan 

Abstract 

Tj\i' 

Jfcf  surveyssome  of  the  more  inportant  theoretical  results  about  the  structure 
of  the  reliability  function  of  a  k  out  of  n  system,  and  indicate^how  these  results 
may  be  used  to  obtain  easily  calculable  bounds  for  the  reliability  of  a  specified 
k  out  of  n  system. 


JU _ Introduct  ion . 


A  system  of  n  components  which  functions  if  and  only  if  at  least  k  of  the 
components  function  is  called  a  k  out  of  n  system.  Parallel  systems  are  1  out 
of  n  systems,  fail  safe  systems  are  n-1  out  of  n  systems,  and  series  systems  are 
n  out  of  n  systems.  We  make  the  assumption  that  the  n  components  in  the  system 
operate  independently  of  one  another.  Therefore  if  £  =  (p^ ,  ....  pn)  is  the  vec¬ 


tor  of  component  reliabilities  and  £  *  (e^,  ...,  en)  represents  any  vector  with 


components  equal  to  zeroes  or  ones. 


then 


\  <R) 


I 


e,  e  2k 

—  l  n 


c\  €t>  1-e 

Pl  -  Pn  (1-PP 


is  the  probability  that  k  or  more  of  the  components  function.  We  refer  to 
,n 


hj^JJ)  :  [0,  1]  -*•  [0,  1]  as  the  reliability  function  of  a  k  out  of  n  system  with 
independent  components.  For  parallel,  failsafe,  and  series  systems,  the  relia¬ 
bility  function  hk(£)  is  a  quite  manageable  function,  even  for  large  values  of  n. 

In  general,  however,  the  behavior  of  the  function  h^O?)  on  [0,  l]n  can  be  quite 
complex,  and  the  calculation  of  {^(£3  at  a  single  vector  of  component  reliabilities 
can  be  quite  cumbersome.  In  the  relatively  simple  case  of  a  5  out  of  8  system  for 
example,  evaluating  h,.(£)  would  usually  involve  calculating  93  products  of  8  num¬ 
bers.  In  this  expository  paper  we  survey  some  of  the  important  theoretical  re¬ 
sults  which  have  been  obtained  about  the  reliability  function  30(1  indicate 

how  these  may  be  used  to  reduce  calculations  and  obtain  good  bounds  for  the 
reliability  of  a  k  out  of  n  system. 


2.  Basic  Theoretical  Results. 

Hoeffding  (1956)  considers  the  problem  of  finding  the  maximum  and  minimum 

n 

of  h.  (p|,  ....  p  )  subject  to  the  constraint  that  £  p.  is  held  fixed.  His  results 

i®l  x 

are  presented  in  terms  of  the  number  of  successes  in  independent  Bernoulli  trials. 
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It  is  clear  however  that  if  S  is  the  number  of  successes  in  n  independent  Bernoulli 
trials  with  respective  success  probabilities  p^,  . ..,  pn<  then  Prob(Sik)  =  h^ 

(Pj,  ...»  Pn) •  For  a  given  vector  £  =  (p^,  ....  pn)  of  component  probabilities. 


_  n  n  n 

we  let  p  =  l  p./n,  [  l  p  ]  be  the  greatest  integer  less  than  or  equal  to  £  p  , 
i=l  1  i=l  1  i-1  1 

n  n 

and  (1,  ....  1,  I  P-  -  (  J  p. ],  0,  ... ,  0)  be  the  vector  with  [£p.]  of  its  n 

«  X  .  .  1  X 

1*1  1=1 

coordinates  equal  to  1.  Hoeffding  proves  in  particular  the  following  theorem: 


Theorem  1.  For  any  vector  £  =  (p^,  pn)  of  component  probabilities  we 


have: 


I  *  h.  (1,  . .  • ,  1,  J  p,  " 


i*l 


n 


I  Pi 
i=l  1 


,0,  ...»  0)  2  h^fp^,  •••»  p^)  ^  h^(p,  . ..,  p) 
n 

if  l  Pi  a  k. 


while 

0  -  hku 


n  I  n  I  _  _ 

,  1,  I  Pi  "  1  Pj  »  °»  •••»  5  hk^pl*  *•**  pn^  *  hk^p* 

i=l  1 i=l  _| 


if  l  P4  5  k-1. 

i=l 


Hoeffding  also  establishes  upper  and  lower  bounds  for  h,  (£}  for  the  case 
n  K 

when  k-1  <  \  p.  <  k;  however  these  bounds  are  considerably  more  complicated. 
i=l  1 

We  present  the  first  of  four  figures  in  an  attempt  to  give  some  geometrical 

insight  into  the  behavior  of  h^(£) .  For  any  a  such  that  0  s  a  <  n,  we  define 

Xa  s  {£:  J  p.  »  a  and  0  <.  p.  £  1  for  isl,...,  n>.  Although  the  convex  set  \a  * 
u  i*l  1  n 

is  the  intersection  of  the  hyperplane  {x:  J  x,  »  a}  with  the  n>dimensional  unit 

i»l 

cube,  it  will  be  represented  geometrically  by  a  "line"  in  the  following  figures. 
The  n-dimensional  cube  itself  is  represented  by  a  "square",  and  by  the  "diagonal" 


we  will  mean  the  line  of  points  all  of  whose  coordinates  are  equal.  Theorem  1  says 
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that  h^(g)  restricted  to  Xo  assumes  a  maximum  of  one  at  the  Extremities '  of 
XQ  and  a  minimum  of  h^Ca/n,  ...»  a/n)  along  the  diagonal  whenever  a  z  k,  while 
it  assumes  a  maximum  of  h^Ca/n,  . ..,  o/n)  on  the  diagonal  and  a  minimum  of  0 


at  the  extremities  of  Xa  whenever  a  s  k-1. 


(Figure  1  goes  here.) 


The  theory  of  majorizstion  and  Schur  functions  provides  an  elegant  format 
for  presenting  many  of  the  results  concerning  the  reliability  function  h^(£) . 

Given  a  vector  x  *  (Xj,  ....  xj ,  let  x^  2  x^  £  . . .  2  x^  denote  a  decreasing 
rearrangement  of  x^ ,  ...,  xn«  (we  write  x_  >m  jr)  if 


2  L  Xm  for  ••••  n-1 


n  n 

Hardy,  Littlewood  and  Polya  (1952)  show  that  x  >m  £  if  and  only  if  there  exists 
a  doubly  stochastic  matrix  n  such  that  £  »  xlt.  If  x  >“  then  the  coordinates 
of  x  are  more  "dispersed'  than  those  of  £  (See  Marshall  and  Olkin  (1979)  for  an 
excellent  treatment  of  majorization.)  Sthur  functions  are  real  valued  functions 
which  are  monotone  with  respect  to  the  partial  ordering  of  majorization.  A 
function  h  with  the  property  that  x  >m  £  *>  h(xj  2  (S)  h(£)  is  called  Schur 
convex  (Schur  concave) .  A  useful  characterization  of  Schur  convexity  (-concavity) 
is  given  by  the  Schur-Ostrowski  condition,  which  states  that  a  differentiable 
permutation  invariant  function  h  defined  on  Rn  is  Schur  convex  (Schur  concave) 
if  and  only  if 

(x4  -  Xj)  (1“^-  |j  )  2  (s)  0  for  all  i,  j  and  xeRn. 

Some  useful  comparisons  for  h^(£)  can  be  obtained  by  working  with  the 
hazard  transform.  The  hazard  transform  z^:  [0,  »)n  [0,  •)  of  h^  is  defined 
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by 

zk(Ri'  •••»  Rn)  =  “1°S  ^lt^e  ^ *  •••*  e  n)  • 

Pledger  and  Proschan  (1971)  prove  the  following  theorem  about  z^: 

Theorem  2.  The  hazard  transform  z^  of  a  k  out  of  n  system  is  increasing 
and  Schur  concave  on  [0,  ®)  n. 

A  particularly  interesting  corollary  of  Theorem  2  is  that  h^p^,  •••»  Pn) 
a  \fPg*  •••»  Pg)  where  pQ  is  the  geometric  mean  ^SjPj)  X  •  Together  with 
Theorem  1  this  implies  that 

^k^P'  •••»  P)  ^  ^(pQ,  ****  P(p 

whenever  n 

[  p.  s  k-1. 

1=1* 

Gleser  (1975)  uses  majorization  to  obtain  refinements  of  the  inequalities 

established  by  Hoeffding.  If  S  is  the  number  of  successes  in  n  independent 

trials  with  respective  success  probabilities  Pj . Pn»  then  Samuels  (1965) 

shows  that  Prob(S=j)  is  unimodal  with  mode  at  [  ?  p.J  or  [  £  p.*l].  Using  this 

i=l  1  i=l  1 

fact  together  with  the  Schur-Ostrowski  condition,  Gleser  proves  the  following 
theorem: 

Theorem  3.  The  reliability  function  hk(£):  [0,  1]°  ■*  [0,  1)  is  Schur  convex 

in  the  region  where  f  p.  z  k*l  and  Schur  concave  in  the  region  where  £  p.  s  k-2. 

i=l  1  i=l  1 

Figure  2  may  help  to  illustrate  Theorem  3.  The  region  shaded  by  vertical 

lines  indicates  where  is  Schur  convex,  while  the  area  shaded  by  horizontal 

lines  indicates  whore  h^(g)  is  Schur  concave. 


(Figure  2  goes  here.) 


If  £  -  (Pj,  pn)  where  n>2,  we  denote  by  £*^  the  vector  in  [0,  l]n  2 
obtained  by  deleting  the  i**1  and  j**1  coordinates  of  £.  For  any  r^O,  we  let 
h^Cj)1*)  denote  the  probability  that  exactly  r  of  the  n-2  components  (with 
respective  probabilities  given  by  g1^)  function.  Boland  and  Proschan  (1983) 
prove  that 

h*_2(Elj)  S  (£)  hk_J(£ii)  for  i*j,  k*2 

whenever 

k-l 

p!l  2  ^0r  a11  Vi.j.  This  resuIf*  coupled  with  the  Schur- 

Ostrowski  condition,  enable  Boland  and  Proschan  to  prove  the  following  Schur 
property  of  h^ (£) : 

Theorem  4.  The  reliability  function  h^(£):  [0,  l]n  -*•  [0,  1]  is  Schur  convex 

in  the  region  [  l]n  and  Schur  concave  in  the  region  [0,  jj^y  ]n. 

k-l  k-l 

The  point  (yy,  ....  yy)  on  the  diagonal  is  a  focal  point  in  distinguishing 

k-l 

the  regions  where  h^(£)  is  Schur  convex  and  Schur  concave.  The  number  y-y  plays 

an  important  role  in  the  study  of  the  function  of  one  variable  h.  (p)  = 
n  n  i  n_j 

I  (i)  P  (1-p)  .  (Note  that  h.  (p)  is  the  restriction  of  ly(£)  to  the  diagonal.) 

^"1  k-l  k-l 

h^(p)  is  convex  on  the  interval  [0,  y-y)  and  concave  on  the  interval  [y-y,l )  (See 

Barlow  and  Proschan  (1965)).  Figure  3  may  help  to  interpret  Theorem  4.  Again 

the  vertically  (horizontally)  shaded  region  indicates  where  hk(£)  is  Schur  convex 

(Schur  concave).  Figure  4  illustrates  the  combined  results  of  Theorem  3  and 

Theorem  4. 


(Figure  3  and  Figure  4  go  here.) 
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3.  Applications. 

We  now  give  a  number  of  examples  in  order  to  demonstrate  how  the  properties 
of  h^CE.)  presented  can  help  us  both  in  comparing  various  k  out  of  n  systems  and 
in  calculating  bounds  for  the  reliability  of  specified  k  out  of  n  systems. 

Example  1.  Let  us  consider  a  3  out  of  4  system.  Theorem  4  in  particular 

implies  that  a  system  with  component  reliabilities  (.7,  .8,  .9,  1.0)  is  superior 

(has  higher  reliability)  than  a  system  with  component  reliabilities  (.75,  .75, 

.95,  .95)  which  in  turn  is  superior  to  one  with  component  reliabilities  (.85,  .85, 

.85,  .85).  Note  that  for  all  three  of  these  systems  the  four  component  reli- 

k-1  2 

abilities  sum  to  3.4,  and  that  each  couyonent  reliability  exceeds  =  -j. 

Theorem  4  also  implies  that  the  system  with  component  reliabilities  (.2,  .3,  .5,  .6) 
is  inferior  to  one  with  component  reliabilities  (.2,  .4,  .4,  .6)  which  in  turn  is 
inferior  to  one  with  component  reliabilities  (.4,  .4,  .4,  .4). 

Example  2.  We  now  consider  evaluating  the  reliability  of  a  5  out  of  8 
system.  As  previously  indicated,  even  for  this  example  where  k  and  n  are  rather 
small,  the  standard  method  for  evaluating  hg(p^,  ...»  pg)  would  normally  involve 
(if  for  example  the  eight  component  probabilities  are  distinct)  calculating  the 
sum  of  93  products  of  8  numbers  each.  We  will  now  indicate  by  judiciously  ’aver¬ 
aging'  some  (or  all)  of  the  component  probabilities,  how  more  easily  calculable 
bounds  may  be  confuted  for  the  reliability  of  the  system. 

Given  a  vector  ■  (p^,  ...,  pg)  of  component  probabilities,  we  lose  no 
generality  in  assuming  that  Pj  5  P2  -  •••  8  Pg*  ^et  118  use  t*'e  following 
notation:  p^  *  (pA  ♦  p^)/2  and  *  (p.  ♦  p^  ♦  p£  ♦  Pm)/4  f°r  any  integers 

i,  j,  H,  m  between  1  and  8.  As  before  p  and  pQ  will  denote  respectively  the 
arithmetic  and  geometric  means  of  the  components  of  the  vector  £  =  (p^ . pg) . 


-  7  - 

Note  that  £  =  (Pj,  P2»  Pj.  P4»  Pg.  P6.  P?»  Pg)  >  (P12»  P12»  P34*  P34»  Ps6*  p56’ 

p78’  p78* 

>  lp1234'  P1234*  p1234*  P1234*  PS678*  P5678*  P5678’  PS678J 

>m  (P»  P.  P.  P,  P.  P,  P,  P)  =  £ 

The  calculation  of  h^(£)  becomes  easier  as  the  number  of  distinct  component  values 

in  £  decreases.  While  the  calculation  of  hg(£)  could  involve  adding  93  products, 

the  calculations  of  h5(p12,  p12,  P34,  P34,  P56<  PS6>  p7g,  P78) ,  hs(Pi234*  P1234’ 

p1234.  p1234*  P5678*  p5678*  p5678*  PS6783  3,1(1  h5 ®  WOUld  involve  addin«  respec¬ 
tively  at  most  31,  10,  and  4  products.  If  the  vector  components  of  £  =  (p  ,  p^, 

k-1  4  4 

p3*  p4’  p5*  p6’  p7*  p8^  are  al1  "  iTT  =  7  fresPectively  *  y)»  then  we  can  easily 
determine  some  lower  (upper)  bounds  for  hg (£)  by  using  Theorems  1  and  4.  To 

illustrate  the  accuracy  of  these  bounds  we  compute  some  specific  probabilities. 

4 

a)  Let  £=  (.60,  .64,  .70,  .74,  .80,  .84,  .90,  .94).  Each  p^  2  y  and  we 
observe  that 


h$(£)  =  .92097 

«=  hg( .60,  .64,  .70,  .74, 

* 

Tf 

00 

o 

oo 

.90, 

.94) 

2  .92069 

=  h$( .62 ,  .62,  .72,  .72, 

.82,  .82, 

.92, 

.92) 

2  .91896 

=  hg(.67,  .67,  .67,  .67, 

.87,  .87, 

.87, 

.87) 

2  .91201 

=  h,. (.77,  .77,  .77,  .77, 

.77,  .77, 

.77, 

.77) 

b)  Let  £  *  (.61,  .62,  .63,  .64,  .65,  .66,  .67,  .68).  Then  again  each 

4 

p^  2  y  and  we  observe  that 

hs(£)  =  .69580  «  h5(.61,  .62,  .63,  .64,  .65,  .66,  .67,  .68) 

2  .69580  =  h,.  (.615 ,  .615,  .635,  .635,  .655,  .655,  .675,  .675) 

2  .69576  ■  hg(.625,  .625,  .625,  .625,  .665,  .665,  .665,  .665) 

2  .69562  •  hg(.645,  .645,  .645,  .645,  ,64S,  .645,  .645,  .64S) 


c)  If  £  *  (-10,  .14,  .20,  .24,  .30,  .34,  .40,  .44),  then  each  component 

4 

p.  ^  y  and  we  observe  that 


=  .03188 

=  h5(. 

10, 

.14, 

.20, 

.24, 

.30, 

.34, 

.40, 

.44) 

.03206 

=  h5c. 

12, 

.12, 

.22, 

.22, 

.32, 

.32, 

.42, 

.42) 

.03318 

=  hs(. 

17, 

•  17, 

.17, 

.17, 

.37, 

.37, 

.37, 

.37) 

.03768 

=  h,(. 

27, 

.27, 

.27, 

.27, 

.27, 

.27, 

.27, 

.27). 

Note  that  for  the  given  vector  £,  the  geometric  mean  of  the  coordinates  is 

Pq  =  .24276.  Using  the  corollary  to  Theorem  2  we  can  compute  the  lower  bound  for 

hk(jj)  which  is  given  by  hgQ)^)  =  .02410. 
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